Collective modes of doped graphene and a standard 2DEG in a strong magnetic field: 
linear magneto-plasmons versus magneto-excitons 
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A doped graphene layer in the integer quantum Hall regime reveals a highly unusual particle-hole 
excitation spectrum, which is calculated from the dynamical polarizability in the random phase 
approximation. We find that the elementary neutral excitations in graphene in a magnetic field are 
unlike those of a standard two-dimensional electron gas (2DEG): in addition to the upper-hybrid 
mode, the particle-hole spectrum is reorganized in linear magneto-plasmons that disperse roughly 
parallel to o) = v^q, instead of the usual horizontal (almost dispersionless) magneto-excitons. These 
modes could be detected in an inelastic light scattering experiment. 
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I. INTRODUCTION 

Particlc-holc excitations in a standard two-dimensional 
electron gas are incoherent and form a continuum in 
the energy-momentum plane This continuum does not 
have a uniform weight but features some structure that 
hints at possible coherent excitations. If one includes 
interactions between the electrons, coherent excitations 
may emerge from the continuum. This is the case of 
the plasmon, which appears once long-range Coulomb 
interactions are taken into account. The plasmon mode 
is long-lived outside the continuum, but it is Landau 
damped once it penetrates it. If instead of including in- 
teractions, one turns on a magnetic field, the electronic 
energy is quenched into Landau levels (LLs), and the 
particle-hole continuum becomes discrete in energy, while 
remaining essentially continuous in the momentum direc- 
tion. When both the magnetic field and Coulomb inter- 
actions are considered, the continuum is reorganized in 
magneto-excitons - transverse excitations that now ac- 
quire a dispersion^ - and an upper-hybrid mode^ exists 
outside what used to be the continuum. 

Graphene, a recently discovered carbon material, is 
attracting a lot of interest due to its unique electronic 
properties (see Ref. [H). Electrons in graphene may 
be viewed as a particular form of the 2DEG. How- 
ever, due to its underlying triangular lattice with a two- 
atom basis, the electrons of graphene are described by 
a massless Dirac equation instead of the usual effective- 
nrass Schrodinger equation. The screening properties of 
graphene are different from that of the 2DEG, as it may 
be seen from the polarization and dielectric functions in 
the two cases i^i^i^ii^ Furthermore, the single-particle 
spectral function reveals the particular chiral proper- 
ties of the Dirac-like quasiparticles , 1 1 ■ Another salient 
physical consequence is a peculiar integer quantum Hall 
effect -i^iii 

In this paper, we show that collective excitations of 
graphene in a strong magnetic field are much unlike that 
of the standard 2DEG [from now on, we use the term 
2DEG to denote a standard 2DEG with a parabolic 



band dispersion, in contrast to graphene]. The mag- 
netic field also reorganizes the particle-hole continuum 
but instead of revealing horizontal lines, diagonal lines 
emerge. When Coulomb interactions are included within 
the random-phase approximation (RPA), the diagonal 
excitations acquire coherence. We shall refer to them 
as linear magneto-plasmons in order to distinguish them 
from the upper-hybrid mode - which is also present and 
which may be considered as the plasmon mode modified 
by the magnetic field - and from the usual (horizontal) 
magneto-excitons, which are blurred in graphene. Sev- 
eral recent theoretical works considered collective excita- 
tions in graphene in a magnetic field, concentrating either 
on magneto-exciton a^^d^d"^ or on the plasmon modcj ^^d^ 
Here, we show that concentrating on horizontal magneto- 
excitons, in spite of its success in the 2DEG^ is not suffi- 
cient to reveal the complete structure of the particle-hole 
excitation spectrum (PHES) in graphene, which is domi- 
nated by linear magneto-plasmons that disperse roughly 
parallel to tj = vpq. 

The paper is organized as follows. In Sec. II we 
calculate the polarization function of a 2DEG and that 
of doped graphene in the integer quantum Hall regime. 
In Sec. HI we discuss the main features of the non- 
interacting PHES of graphene and a 2DEG, and the 
different nature of the collective modes in each case, 
once electron-electron interactions are taken into ac- 
count. The conclusions are summarized in Sec. IV. 



II. POLARIZATION FUNCTION 

The polarization operator H(q, cj) may be viewed as 
the particle-hole pair propagator, the poles of which yield 
the dispersion relation and damping of the collective ex- 
citations. Its imaginary part is related to the dynamical 
structure factor by 
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which plays the role of a spectral function for the particle- 
hole propagator. The PHES is defined as the (w, q) region 
of non-zero spectral weight S{q,uj) ^ 0. Peaks in the 
spectral density are interpreted as collective excitations; 
their dispersion relation and damping may be extracted 
from the position and the width of the peaks (see e. g. 
§3.2.7. in Ref. Q). 

The Fourier-transformed polarizability for non- 
interacting electrons (we use a unit system with fi, = 1), 

n°(q, u;)=J 1^ Tr[GO(k, u;')G'ik + q, ^' + c.)] , (1) 

is given in terms of the free single-particle Green's func- 
tions G°(q, a;). Whereas, as we review below, the single- 
particle Green's functions for the conventional 2DEG are 
simple scalar functions, those of graphene are 2x2 ma- 
trices due to the two sublattices A and B. The trace (Tr) 
in Eq. ([1]) is, thus, relevant only in the case of graphene. 



A. Polarizability in the 2DEG 

The Hamiltonian for electrons in the conventional 
2DEG is simply that of free electrons with a band mass 
rrib, Tl = 7r^/2mb, where the gauge-invariant momentum 
operators, tt ~ p + eA(r) takes into account the coupling 
to the magnetic field B = V x A, where A is the vector 
potential. The operators tt may be expressed in terms of 
the usual ladder operators with the help of 



which satisfy the commutation relations [a, a^] — 1, and 
Ib = l/^feB is the magnetic length. One obtains the 
usual eigenstates, |n, m), and the eigenvalue equation 
7i|n,m) = ujc{n + l/2)|n,TO), where uj^ ~ l/miylg = 
eB I rub is the cyclotron frequency. The result is the usual 
LL spectrum where the levels are labeled by the quan- 
tum number n, the eigenvalue of a^a. The additional 
quantum number m is associated with the guiding cen- 
ter operator, and varies between and Nb — 1, where 
Nb = A/2TTlg, in terms of the sample area A. 

With the help of the eigenstates \n,m), one may ex- 
press the Green's functions G°(q, w) for non- interacting 
electrons in Fourier space, as 

(q|n,77i)(n,77i|r = 0) 
r::^, +J(5sgn(^„) 
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where ^„ = ujc{n + 1/2) — ep is the energy of the n-th 
LL with respect to the Fermi energy ep, which we choose 
to lie between two LLs of the conduction band (inte- 
ger quantum Hall regime), and S accounts for disorder- 
induced level broadening {6 0+ in the clean limit). 

Eq. ([2]) allows us to calculate the polarization function 
(H]) for the 2DEG, 
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where Np is the index of the last occupied LL, which is 
fixed by the filling factor, and uj^ — * a;~ indicates the 
replacement uj -\- i5 ^ —lo — i6. The form factors for the 
2DEG readi 
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with n> = max{n, n'} and n< = min{ri, n'}, and L™ are 
associated Laguerre polynomials. 

B. Polarizability in graphene 

For graphene in a magnetic field, the electronic Hamil- 
tonian may be written as^ 



H = 



V2' 



a 
at 



(4) 



where v-p = itacc/^ is the Fermi velocity, in terms of the 
the nearest-neighbor hopping integral t ~ 3eV, and the 
carbon-carbon distance Occ — 1.4A. Strictly speaking, 
Eq. (HI) is only valid on one of the valleys, namely if, 
that of the valley K' being — 7i*. However, we concen- 
trate, here, only on processes that do not couple the two 
different valleys, such that a discussion of the Hamilto- 
nian ^ is sufficient, and the twofold valley degeneracy 
may be accounted for by a simple factor of 2. 

The Hamiltonian ^ yields the relativistic LL spec- 
trum 



Ae„ = A — A 

Lb 
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where A denotes the band index, A = -f for the conduc- 
tion and A = — for the valence band. The associated 
eigenstates are the 2-spinors 
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\n ~ 0, m) 



for n> I 
for n = 



where m = 0,1, Nb — 1, and \n,m) are the corre- 
sponding eigenstates of the Hamiltonian with quadratic 
dispersion, introduced in the previous subsection. Due to 
the 2-spinor structure of the wavefunctions in graphene, 
the associated single-particle Green's functions are 2x2 
matrices, 
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with a = A{B) for electrons on the A{B) sublattice, 
C = -|-(— ) for electrons in the K{K')-va.\lcy, and k is 
the electron momentum measured from the Dirac points, 
au,.. In the denominator, we have defined 
A^„ = Ae„ — ep, the energy difference between the level 
Ae„ and the Fermi energy ep. Furthermore, we neglect 
Zeeman and possible valley splittings. The matrix /An(q) 
for the if-valley is 
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/K;An(q) 



Ef l*2(q|n - l,TO)(n - l,TO|r = 0) 
UAi;2;(q|n,m)(n-l,m|r = 0) 



-iAl*2* 
2f ( 



(q|n — 1, m){n, m\r = 0) 
q|n, to) (n, TO|r = 0) 
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where we have introduced the simphfied notation 1* = 
y/il - (5„,o)/2, 2* = V(l + '5„,o)/2. The bare polariza- 
tion bubble is 



Nf 



n«(q, u;) = ^ (q, u;) + n™^(q, uj) , 



(7) 



n=l 



where the index Np of the last occupied LL is now re- 
lated to the filling factor ly hy v = ANp + 2 due to the 
fourfold spin and valley degeneracy. There is a non-zero 
vacuum polarization, n'"'^(q, w) = — X]^=i n^""^(q, w), 
the vacuum corresponding to undoped graphcne. Here, 
Nc is an ultra-violet cutofS^ and 



n-l 



In the previous expressions, we have used 
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- + iuj+^u;-), (9) 



A' n'=0 

-I- ^ ^ n^^/,(q, + n^~^(q,cj). (8) where the form factor is 



A' n'—7i-\-l 



Notice that Eq. (fTU]) agrees with the form factor obtained 
in Ref. but not with those obtained in Rcf. ilSiiT^. 
Because of the presence of nodes in the LL wavefunctions 

(zeros of the Laguerre polynomials), J^„„/(q) contains 
zeros and, therefore, so docs ImH^ as a function of q at 
fixed Lu. These zeros will play an important role when 
discussing the structure of the PHES in Sec. III.C. Note 
that n^(q, uj) = — n^^(q, uj) and that the n = LL does 
not contribute to the polarization. 



in graphene has linear edges and contains two different 
regions because of intra- (I) and inter-band (II) processes 
[see Fig. l(a)| . The spectral weight in graphene is con- 
centrated around the diagonal uj = vpq due to the chiral- 
ity factor [1-1- A A' cos 6']/2fii which is the B = equivalent 
ofEq. pO]) . 



Particle-hole excitation spectrum for 
non-interacting electrons 



III. RESULTS 



The dynamical polarization function of graphene, 
^^raph^ is 1^0^ compared to that of the 2DEG, ^2DEG'^ 
For B = 0, ^'^raph been calculated first in the context 
of intercalated graphite^ and later for doped graphene^^ 
resulting in a particle-hole continuum different from that 
of a 2DEG. The edges of the zero-field PHES for graphene 
and a 2DEG are drawn as black lines in Fig. [TJ While 
the PHES for a 2DEG is made of a single region limited 
by linear-parabolic boundaries [see Fig. [l](b)], the PHES 



In the presence of a magnetic field, the largest con- 
tribution to the polarization also comes from the re- 
gion around w = vpq [see Fig. [1] (a)]. In addition, 
ImH°^^p^(q, cj) is finite not only close to the diagonal. 



but also along regions above and below lu 



diagonal yellow stripes in Fig. 

Fig. 1(b) we show the corresponding density plo^ of 



vpq [see the 
For comparison, in 



Im H2]3gQ(q, uj) calculated for the same Fermi momen- 
tum kp ^ V^Np + I/Ib- 

These results reveal the main differences between the 
PHES of graphene and that of a 2DEG in the presence of 
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(c) (d) 

FIG. 1: (Color online) Density plot of the spectral function Imn(q,a;) for doping Np = 3, ultraviolet cutoff Nc = 70 and disorder 
broadening <5 = 0.2vp/lg for grapliene [panels (a) and (c)] and 5 = 0.2aJc for a standard 2DEG [(b) and (d)]: in the absence of interactions 
[(a) and (b)] and including interactions in the RPA [(c) with rs ~ 1 and (d) with rs ~ 3]. The solid lines indicate the zero field limits of 
the PHES and arrows point to the islands discussed in the text. In panel (a), label I (II) indicates the intra- (inter-) band region of the 
PHES in graphene. 



a magnetic field. While the PHES for a 2DEG [Fig. [^b)] 
features a set of well-defined horizontal lines, with a slight 
modulation parallel to the boundaries of the particle-hole 
continuum, that of graphene is remarkably different. In- 
deed, its PHES [Fig. [Ija)] is dominated by this modula- 
tion, visible as bright regions in the diagonal, whereas the 
horizontal lines are now hardly visible, as a consequence 
of the non-equidistant LL spacing and the presence of a 
finite disorder broadening (5, as discussed in more detail 
below [see Sec. IIII Cj . The PHES is, thus, discretized into 
diagonal lines (almost parallel to uj = v-pq) in the (q, lo)- 
plane. For Fig. [1] we have chosen S ~ 0.2vf /Ib, which 
is a realistic value for the currently realized graphene 
samplesi^ Lower values of 6 lead to a clearer definition 
of the horizontal lines [see Fig 



2(a) 



low energy region of the PHES in a c 



for a zoom of the 
eaner system]. 



B. Particle- hole excitation spectrum for interacting 
electrons (RPA) 

Electron-electron interactions yield coherent modes 
that emerge, both in the 2DEG and in graphene, from 
the regions of highest weight in the PHES. Within the 
RPA, the renormalized polarization function is given by 



H 



RPA 



n"(q,c.) 



l-i;(q)HO(q,o;)' 



where w(q) = 27re^/£b|q| is the unscreened 2D Coulomb 
potential, Sb the dielectric constant, and inter- valley pro- 
cesses, which are suppressed in u/Ib, are neglected i2i 
The results are shown in Fig. 1(c) 2(c)| for graph ene 
(with Ts ' ' 



a 2DEG (with 



e'^/ebVF ~ 1) and in Fig. |l(d)j |2(d)| for 



2mbe / ebk[ 



1 or 3). In the 



2DEG, Coulomb interactions lead to the appearance of 
dispersive magneto-excitons [see Fig. |2(d)| , whereas in 
graphene, the diagonal lines of the non-interacting PHES 
[Fig. HJa)] become coherent collective modes: linear 
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(c) (d) 



FIG. 2: (Color online) Density plot of Imn(q,a;) for the low energy region, for doping Np = 3 and ultraviolet cutoff Nc = 70. Panels (a) 
and (c) correspond to graphene, whereas (b) and (d) correspond to a 2DEG. The disorder broadening S = 0.05 in units of vp/lg and lUc, 
respectively, and the interaction strength ~ 1 in panels (c) and (d). Notice the different energy scale with respect to Fig. [T] Arrows 
indicate islands as discussed in the text. 



magneto-plasmons that arc now clearly visible as peaks 
in the spectral function. They disperse roughly parallel 
to (jj = vpq and are more pronounced in the inter-band 
region of the PHES [see Fig. [T^c)]. Notice that even in 
cleaner samples [S = 0.05vf/Ib, Fig. 2(c) , the horizon- 
tal structure due to the LL quantization is highly sup- 
pressed in Imn^^'^(q, cj). Furthermore, there has been 
a transfer of spectral weight from the long-wavelcngth 
region of the non-interacting PHES to the upper-hybrid 
mode that starts dispersing in the gapped region of the 
spectrum [see Fig. 2(c) . In a 2DEG, this upper- hybrid 
mode^ can be seen as the plasmon mode modified by 
the magnetic field. It is a plasmon-cyclotron collective 
mode that has a dispersion relation uj = [uSl + tj^]^/^, 

where Wp ~ ^ epe^g/eb in the long-wavelenght limit. As 
it is the most intense mode, Its square-root behavior is 
clearly visible in Fig. [Hd) [in Fig. 2(d) it only appears as 
a maximum in the magneto-cxcitons dispersion relation] 
and also in Fig. 1(c) 2(c) and 3(b) (d) for graphene. 



C. Discussion 

The structure of the PHES in graphene and the 2DEG 
may be understood in the following manner. First, one 
notices that the boundaries [black lines in Fig. [1] are 
due to the momentum conservation of the electron-hole 
pairs. The boundaries reflect the different B = be- 
havior of the dispersion relations in graphene (linear in 
q) as compared to the 2DEG (quadratic in q). Second, 
for _B 7^ 0, the zero-field continuum becomes chopped 
into horizontal lines. The main difference between the 
PHES in the 2DEG and that in graphene stems from the 
different LL quantization. 

In the case of the 2DEG, one notices a constant spac- 
ing of the horizontal lines, lo ~ mcuc- There are several 
contributions to each horizontal line from the allowed LL 
transitions, n n' , with m ^ n' — n. For a given inte- 
ger m, the number of contributing LL transitions is the 
minimal value minjiVi? -t- 1,to} of Np -\- 1 and m. The 
lowest-energy horizontal line (m ~ 1), thus, arises from 
a single LL transition. Np Np + 1. In this case, the 
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FIG. 3: (Color online) Density plot of Im n^^'*(q, a;) for graphene, for disorder broadening 5 = 0.2vp/lg [plots (a) and (b)] and 
S = 0.05vp/lg [plots (c) and (d)], and doping Np = 1 in plots (a) and (c), and Np = 6 in plots (b) and (d). 



line consists of Np + 1 well-separated islands [see arrows 
in Figs. l(b)| and |2(b)| , which is the name we give to 
regions of high spectral weight, and which are separated 
by regions of low spectral weight reflecting the nodes of 
the LL wave-functions. For larger values of to, several 
LL transitions n — > n' contribute to the same horizontal 
line. Each of these transitions has n -I- 1 islands. Because 
this number is different for each transition, the islands 
that arise from different transitions overlap and fill in 
the region of low spectral weight. As a result, the hori- 
zontal lines appear continuous [compare the to = 1 and 
3 horizontal line in Figs. 



to 



1(b) and 2(b) 



and the 

node structure of the LL wave-function is visible only in 
form of faint shadows parallel to the boundaries of the 
PHES [see Fig. [1(b)] and [1(d)] . 

In graphene, the relativistic LL quantization gives rise 
to horizontal lines at energies 



(11) 



XV2n] 



where we have assumed that the Fermi energy lies in the 
conduction band, an therefore A' = 1. In contrast to the 



2DEG, with equally spaced LLs, only a single transition 
A, n — > n' contributes to each horizontal line. As a conse- 
quence, the node structure of the LL wavefunctions, and 
therefore that of the form factor ([T0|) . is clearly visible 
in every horizontal line and not only in the one of low- 
est energy [see Fig. 



2(a) 



This explains the presence of 
n + 1 islands in the horizontal line corresponding to the 
71 — > n' transition. These nodes give rise to well-defined 
diagonal dark lines of low spectral weight parallel to the 
boundaries of the particle-hole continuum. Another con- 
sequence of the relativistic LL structure is that the den- 
sity of horizontal lines increases with energy, contrary to 
the case of the 2DEG, where all lines are separated by th e 
constant cyclotron energy [compare Figs. 2(a) and [2(b)] . 
This clear separation of islands at any energy, together 
with the rather large number of horizontal lines and their 
finite disorder-induced width, leads to a stacking of the 
islands of different energies. This is the origin of the ap- 
pearance of diagonal regions of strong spectral weight [see 
Figs 1(a) and 2(a) . The overlap of islands of different en- 



ergies is stronger for higher values of disorder in the sam- 
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pie (increased LL broadening 6) as well as for higher fill- 
ing factors. The relative energy separation between the 
horizontal lines will be smaller for larger Np, because in 
the single-particle graphene spectrum, the density of LLs 
increases with energy. Electron-electron interactions turn 
these original regions of strong spectral weight into co- 
herent excitations (the upper-hybrid mode and the linear 
magneto-plasmons). In fact, we have studied the PHES 
of graphene for different values of doping (i. e. Np 7^ 3) 
and disorder broadening S and found that the physical 
picture exposed above is unaltered, [see Fig. [3]. As a 
general rule, disorder and doping favor the emergence 
of linear magneto-plasmons in graphene as compared to 
the horizontal magneto-excitions. For the experimentally 
relevant values of S and Np, linear magneto-plasmons 
therefore dominate the PHES of graphene in the inte- 
ger quantum Hall regime. It is worth noticing that, at 
fixed carrier density, increasing the filling factor Np) 
is equivalent to effectively lowering the magnetic field. 
Although Fig. [3] represents the PHES at fixed values of 
vp/Ib, i- e. fixed magnetic field, one easily sees that 
at larger values of Np, most of the spectral wei ght i s 
concentrated in the upper- hybrid mode [see Fig. |3(d)| . 
This mode is the only collective excitation in the integer 
quantum Hall regime that evolves continuously into the 
plasmon mode at zero field. In contrast to the upper- 
hybrid mode, the linear magneto-plasmons become less 
intense in the large- A''^? limit and evolve for B ^ {) into 
the incoherent particle-hole continuum. 

IV. CONCLUSIONS 

In conclusion, the magnetic-field particle-hole excita- 
tion spectrum in graphene has been investigated, and 



the results have been compared to those of a standard 
2DEG with a parabolic band dispersion. Most saliently, 
the particular LL quantization in graphene yields lin- 
ear magneto-plasmon modes, which are not captured 
in the usual magneto-exciton approximation. This is 
due to the fact that a single particle-hole process con- 
tributes to a given energy in graphene, whereas for a 
2DEG, there are, in general, many processes contribut- 
ing to the same energy. As a consequence, the node 
structure of the LL wavefunctions in graphene is evi- 
dent at any energy of the PHES, leading to a partic- 
ular spectrum where the relevant modes are diagonal, 
dispersing roughly parallel to v^q. The dominant role 
of linear magneto-plasmons is enhanced by doping and 
Landau level broadening due to disorder. In addition to 
linear magneto-plasmons, electron-electron interactions 
give rise to the upper-hybrid mode as well, which is the 
plasmon mode rcnormalizcd by the magnetic field. This 
alternative scenario of linear magneto-plasmons, as op- 
posed to horizontal magneto-excitons, may find an exper- 
imental proof in the framework of microwave absorption 
or inelastic (Raman) light scattering, which has proven 
to be a powerful tool for measuring the collective excita- 
tions of the 2DEG in a strong magnetic field, see e. g. 
Ref. [13. 



Acknowledgments 

We thank Y. Gallais, P. Lederer, and F. Picchon for 
useful discussions and acknowledge funding from the 
ANR under Grant No. ANR-06-NANO-019-03 and "Tri- 
angle de la Physique" . 



^ G. F. Giuliani and G. Vignale, Quatum Theory of the Elec- 
tron Liquid (CUP, Cambridge, 2005). 

^ D. Pines and P. Nozieres, The Theory of Quantum Liquids 
(W. A. Benjamin, Inc., New York, 1966). 

^ G. Kallin and B. I. Halperin, Phys. Rev. B 30, 5655 (1984). 

* K. W. Chiu and J. J. Quinn, Phys. Rev. B 9, 4724 (1974), 
T. H. Styx, Theory of plasma waves, Mc Graw-Hill, 1962. 

^ A. H. C. Neto, F. Guinea, N. M. R. Peres, K. Novoselov, 
and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 
K. W. K. Shung, Phys. Rev. B 34, 979 (1986). 
J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Nucl. 
Phys. B 424, 595 (1994). 

* B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New Jour- 
nal of Physics 8, 318 (2006). 

® E. H. Hwang and S. D. Sarma, Phys. Rev. B 75, 205418 

(2007) . 

^° J. Sabio, J. Nilsson, and A. H. C. Neto, Phys. Rev. B 78, 
075410 (2008). 

M. Polini, R. Asgari, G. Borghi, Y. Barlas, T. Pereg- 
Barnea, and A. H. MacDonald, Phys. Rev. B 77, 081411 

(2008) . 



^2 E. H. Hwang and S. D. Sarma, Phys. Rev. B 77, 081412 
(2008). 

K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature (London) 438, 197 (2005). 
Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 
(London) 438, 201 (2005). 

A. Iyengar, J. Wang, H. A. Fertig, and L. Brey, Phys. Rev. 
B 75, 125430 (2007). 
" Y. A. Bychkov and G. Martinez, Phys. Rev. B 77, 125417 
(2008). 

K. Shizuya, Phys. Rev. B 75, 245417 (2007). 

M. Tahir and K. Sabeeh, J. Phys.: Gondens. Matter 20, 

425202 (2008). 

O. L. Berman, G. Gumbs, and Y. E. Lozovik, Phys. Rev. 
B 78, 085401 (2008). 
2° T. Ando, J. Phys. Soc. Jpn. 76, 024712 (2007). 

M. O. Goerbig, R. Moessner, and B. Dougot, Phys. Rev. 
B 74, 161407(R) (2006). 

M. A. Eriksson, A. Pinczuk, B. S. Dennis, S. H. Simon, 
L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 82, 2163 



8 



(1999). 

The continuum approximation is valid up to an energy 
~ t. This leads to a cutoff in the Landau level index 
Nc ^ 10*/-B[r]. As the LL separation in graphene de- 
creases with n, it is possible to have semi- quantitative re- 



sults from smaller values of N^. Here, for numerical calcu- 
lations, we typically take Nc « 70. There are no noticeable 
differences for larger values of Nc- 



